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Abstract 

We extract the S'-matrix element for dipole-proton scattering using the data on diffrac- 
tive electroproduction of vector mesons at HERA. By considering the full t dependence 
of this process we are able to reliably unfold the profile of the S'-matrix for impact pa- 
rameter values h > 0.3 fm . We show that the results depend only weakly on the choice 
of the form for the vector meson wave function. We relate this result to the discussion 
about possible saturation effects at HERA. 

Introduction 

One of the key issues in deeply inelastic scattering and deeply inelastic diffraction is whether 
parton saturation effects are present in the HERA energy regime. The answer to this question 
has important implications for our understanding of the very early stages of relativistic heavy 
ion collisions where parton saturation would result in the production of a high density and 
high field strength, F^,j ~ l/^fal , state of QCD. Such a state would be a new and exciting 
regime of nonperturbative but small coupling QCD. 

Parton saturation can be viewed as a manifestation of unitarity limits being reached. 
However, unitarity limits are difficult to see directly in deep inelastic scattering and diffrac- 
tion. Thus, most studies of parton saturation at HERA have been done in terms of models 
which explicitly impose unitarity at moderate and small x. The most widely used such 
model is that of Golec-Biernat and Wiisthoff Q where the S-matrix for the scattering of a 
dipole of separation r on a proton is given hy S = e~^ Qs/'^ with the saturation momentum, 
Qg, depending on the energy of the scattering. The fact that the Golec-Biernat Wiisthoff 
model and its generalization work so well at HERA is perhaps some evidence that deep in- 
elastic scattering, and especially diffraction, have reached unitarity limits. 



It would be nice to see the attainment of unitarity limits directly without having to use 
a model. This is not an easy task. Recall that in proton-proton scattering the energy depen- 
dence of total and elastic cross sections changes little between ISR and Tevatron energies. 
The near saturation of unitarity limits was found in proton-proton scattering only when 
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0194, by the Polish Committee for Scientific Research grants Nos. KBN 2P03B 120 19, 2P03B 051 19, 5P03B 
144 20 and by the US Department of Energy. 
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Amaldi and Schubert Q evaluated the proton-proton elastic scattering amplitude in impact 
parameter space and found near blackness for small impact parameters. 

Vector meson production seems to be the best process from which to extract the dipole 
proton elastic scattering amplitude. For example, it is easy to check that our procedure (see 
Eqs. (^)-(|T^) below) would not work for the total diffractive cross section because there is 
not a single definite state produced, a condition which seems necessary in the way we proceed. 

In this paper we copy the Amaldi-Schubert analysis in order to determine the dipole- 
proton scattering amplitude in impact parameter space in terms of diffractive p— production 
data at HERA. Recall that diffractive electroproduction of a /)— meson can be viewed in terms 
of the following sequence of transitions: (i) The virtual photon goes into a quark-antiquark 
pair (dipole). (ii) The dipole scatters elastically on the proton, (iii) The dipole then goes 
into a /9— meson. Thus /o— electroproduction naturally carries information on the dipole-quark 
elastic scattering amplitude. The main uncertainty is in the wavefunction of the /)— meson 
which appears to be reasonably well constrained by previous phenomenology. 

The dipole-proton S'-matrix which we extract as a function of impact parameter is av- 
eraged over the dipole sizes which naturally appear in the 7* and p— meson wavefunctions. 
This averaging does not affect the fact that for weak interactions S is near 1 and for strong 
interactions S is near zero. Thus the smallness of S is a measure of how close one is to 
the unitarity limit, S = 0. Perhaps the best way to measure how close the dipole-proton 
scattering is to the unitarity limit is to note that (1 — S*^) > 1 — (S)^ gives the probability 
that a dipole passing the proton will induce an inelastic reaction at the impact parameter 
in question ((• • •) denotes the average over dipole sizes as explained in Eq. (0)). From our 
analysis this probability is likely significantly greater than 1/2 for < 2GeV^ and for b 
near zero (see Fig. although better large momentum transfer data would be needed to 
definitively determine S* at 6 ~ 0. Thus, it is likely that saturation effects are important in 
this low regime. Our estimate of is about 1 — 1.5 GeV^ ai b ^ 0.3 fm, although this 
estimate has large uncertainties because it depends on knowing the average dipole size as a 
function of Q^. Finally, we determine the total dipole-proton cross sections for x ~ 10~^ to 
be 14.9, 10.6 and 7.5 mb at = 0.45,3.5 and 7GeV^ respectively. 

The outline of the paper is as follows: in section |l| we establish the relation between the 
5-matrix for dipole-proton scattering averaged over the distribution of dipole sizes and the 
diffractive differential cross section for electroproduction of vector mesons. Due to the partic- 
ular properties of this distribution, this formula translates into a method of determining the 
S'-matrix profile in impact parameter space for a dipole of given size r. The only theoretical 
input needed is the vector meson wave function. In section || we discuss the parametrizations 
for the wave functions that we used, and we study in detail the model independence of our 
procedure. In section ^ we apply this method to the analysis of the available HERA data 
and present our main results on the S-matrix profile in impact parameter space. We discuss 
the uncertainties due to the experimental errors on the measured cross section and due to 
the lack of data in the high t region. We also discuss the possible implications for studies of 
saturation effects at high energies. Finally, in section ^ we state the summary of our results. 
We append some more technical points. 
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1 Dipole-proton 6'-matrix and diffractive meson production cross section 



We consider the process of diffractive production of a vector meson, j*—p V—p, shown in 
Fig. 1^. We adopt here the dipole formulation ^ in which the virtual photon 7* fluctuates 
into a qq pair (dipole) which subsequently interacts elastically with a proton, and eventually 
forms a meson bound state. This picture is justified at high center-of-mass energy of the 
-f*p system, since in this regime it is guaranteed that these successive processes are clearly 
separated in time. We assume that the energy is sufficiently high so that s-channel helicity 
conservation holds to a good accuracy, and we limit ourselves to the case of longitudinal 
vector meson production. The kinematics of this process is shown in Fig. ||: we define the 
2-momentum transfer A related to the Mandelstam variable t hy t = — A^. = —q^ and 
X = Q^/ {2P-q) are as usual the virtuality of the photon and the Bjorken variable, respectively. 

In this framework, the amplitude Aq\{x, A,Q) for the scattering can be written in the 
following factorized form: 

A,Q) = d\dz^^^t{z,v-Q)Af-^{x,v,A)^\;^{z,v) , (1) 

h,h 

where A'^~^{x, r. A) is the elementary amplitude for the elastic scattering of a dipole of size 

r on the proton. V'^l^ is the photon light-cone wave function projected onto a state made of a 
quark-antiquark pair of charges Cq, masses niq and respective helicities h and h. This function 
is computed in first order light-cone perturbation theory of quantum electrodynamics, and 
for the longitudinally polarized virtual photon, it reads [^, ^ 

^'/{z,r-Q) = 6,^_-, -|= :^2z{l-z)QKo{er) , (2) 



where 

= Q^z{l- z)+ml . (3) 

Similarly, ^l^y^ is the wave function of the vector meson produced in the final state. We 
discuss it in more detail in the next section. The variable z in Eq. (|T|) is the fraction of 
longitudinal momentum of the virtual photon carried by the quark. 

The amplitude is normalized in such a way that the differential cross section for the full 
process is 

^ = ^K,(x,A,e)r (4) 

Furthermore, the elementary amplitude A'^'y ^(x,r. A) can be related to the S-matrix 
element S{x, r, b) for the scattering of a dipole of size r at impact parameter b |^] 

Alf-P{x, r. A) = J d^h Alf-P{x, r, bje**'^ = 2 j d^h[l- r, b)]e^^^ . (5) 

We can check briefly the consistency of this formula, which defines S'(x,r, b) (see Ref. |^). 
In the standard normalizations adopted here for the amplitudes, the optical theorem takes 
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the form o"^^^ ^(x, r) = ZmiA'^'^ ^(x, r, A = 0). Using formula (|5|), this relation translates into 
the following expression for the total dipole-proton cross section: 

aZ'^'ix, r) = 2 y d^h (1 - Tie S{x, r, b)) . (6) 

On the other hand, taking the amplitude @ squared, dividing it by the flux factor and 
integrating it over phase space, it is easy to see that the elastic cross section is given by 

<f-^(x,r) = I d'h \1-S{x,v,h)\' . (7) 

When S = 0, the unitarity limit is reached which is equivalent to the scattering on a black 
body. One sees from formulae (^) and (0) that the elastic cross section is half the total one, 
as it should be. This justifies the consistency of the normalization adopted for S{x, r, b) in 
Eq. (|5|). We will assume in the following that iA"^ ^(x,r. A) is purely imaginary, i.e. that 
S{x, r, b) is real. 

The main goal of this paper is to extract S{x, r, b) taking advantage of the present exper- 
imental knowledge of the cross section for diffractive vector meson production, see Eq. (||). 
To this aim we express the amplitude which appears in Eqs. (0,^) by means of S{x, r, b) using 
relation Then we take the inverse Fourier transform of the square root of Eq.@, which 
gives 



/ ^VW'"^^'' = E / d'rdzi;';'}*{z,r,Q)2[l-S{x,rM^'v\z,r) . (8) 

h,h 



In the following we suppress the angular dependence of the S'-matrix since in this process one 
is only sensitive to the quantities which are angular averaged. We single out the (logarithmic) 
distribution of dipole sizes at the photon vertex, which is the overlap between the photon 
and the meson wave functions 

p(r,Q) = 27rr2 ^ r; Q) r) , (9) 

h,h 

and we denote N{Q) its integral 

f'°° dr 

N{Q)^ -p{r,Q). (10) 
Jo r 

We then define the mean value of a given function /(r) with respect to the probability measure 
p{r,Q)/N{Q) by 

We now rewrite Eq. (^) using these definitions (^), (p^), (pTj), and we obtain the following 
formula: 
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which shows that the measurement of the differential cross section da/dt enables us to de- 
termine the S-matrix at fixed impact parameter b, averaged over the dipole size distribution 
p{r,Q) defined by Eq. (^). 

One should stress that N{Q) which appears in Eq. (|12|) is the only source of theoret- 
ical uncertainty in this formula for the average S'-matrix. This quantity depends on the 
parametrization of the wave function of the vector meson. However, as we shall discuss in 
detail in Sec. N{Q) is very well constrained and is hardly model-dependent. We now have 
to investigate the meaning of this average over dipole sizes, by studying p{r, Q) in more detail. 

The distribution p(r, Q) has some interesting properties which are quite general and rather 
independent of the particular model for the meson wave function: it is sharply peaked at a 

specific value of r, which depends on like rq = Aj \J Q"^ + "ly and its width, roughly 

independent of Q^, is of order unity on a logarithmic scale. To illustrate these properties, 
we choose a particular model for il)y{z,r) (details will be given in the next section), and 
we plot p(r, Q) for various values of in Fig. |2|. One can easily obtain a rough estimate 
of rq considering the fact that on one hand, the photon wave function behaves like logr at 
small r, and on the other hand, it can be approximated by exp(— er) in the asymptotic region 
r — > oo, see Eq. (||). Assuming furthermore that the vector meson wave function is smooth 
in r, the integrand p(r, Q)/r ~ rKQ{£r) is then peaked around the maximum of the function 
rexp(— er), i.e. around rq ~ 1/e. In the case of longitudinal vector meson production, the 
dominant contribution comes from symmetric configurations for which z is close to 1/2. This 
leads to the above-quoted formula for rq, with A ~ 2. 

When {S{x,r,b))p is significantly below 1, the quantity rq can be interpreted as the 
mean value of the sizes r of the dipoles which participate in the interaction. In this case, the 
dipole-proton amplitude A'^'y~^{x,r,b) is large for any dipole of size r present in the initial 
wave function and thus does not filter these dipoles: the dominant dipole sizes involved in 
the interaction are then effectively distributed according to p(r, Q). However, in the opposite 
case when {S{x,r,b))p ~ 1, the amplitude A'^^^~^{x,r,b) is small for the relevant values of r. 
In this situation, we are sensitive to the behaviour of the dipole cross section at small r, i.e. 
A'^~^{x,r,b) ~ r^. This is known as the colour-transparency property. Then the interacting 
dipoles are effectively distributed according to r^p(r, Q). This results in the dominant dipole 
sizes being shifted towards larger values^ of r ~ for which an estimate can be provided 

using the same technique as previously, and yields ~ 6/ Q"^ + my. 

According to this discussion and considering the fact that one is a priori more interested 
in the case when {S{x,r,b)) < 1, one can make the following approximation: 

{S{x,r,b))p c:i S{x,rq,b) with rq ~ , = , ^~2. (13) 

Thus for fixed and energy Eqs. ( p!^Jl^ ) enable us to determine S{x,rq,b), which is the 
S'-matrix element for the scattering of a dipole of size rq on a proton, at energy Q'^/x and 
impact parameter b. A novel feature of our analysis is the fact that we are able to extract the 
profile of the S'-matrix (or equivalently of the dipole-proton amplitude) in impact parameter 
space. The previous studies of Refs. Q, |, 0, |ll|, p], |l|, |l|, [l|] provided an estimate of the 
dipole cross section integrated over the impact parameter. 

^The quantity Va in this case is the so called "scanning radius" introduced in [^. 
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2 From mesons to dipoles 



As already mentioned earlier in Sec.l the only theoretical uncertainty in the extraction of 
^-matrix profile in impact parameter space is the form of the vector meson wave function 
which occurs in the overlap function p{r,Q) (Eq. (^)) and in its integral N{Q) (Eq. (|l2|)). 
In this section, we present the different models that we choose for the vector meson wave 
function, then we compare them phenomenologically, and finally we evaluate the uncertainty 
on S{x, rg, b) induced by the freedom of choice of the vector meson wave function. 



2.1 Models for the vector meson wave function 



Several different models for the vector meson wave function exist in the literature [H, 0, IC 



16, 17]. All these approaches use the information from spectroscopic models of long distance 
physics. In these models one assumes that the meson is composed of a constituent quark and 
antiquark which move in an harmonic oscillator potential. This results in a wave function 
which has a gaussian dependence on the spatial separation between the quarks. Additionally 
these models are supplemented by the short-distance physics driven by the QCD exchange 
of hard gluons between the valence quarks of the vector meson. Finally, a relativization 
procedure has to be applied. 

Of course, there are many uncertainties in the above procedure of obtaining the wave 
functions for the vector mesons, however they are constrained by model independent features. 
First of all, ipy^ has to satisfy the following normalization condition |Q 

1 = ^ / d\dz\^'/(z,r)f . (14) 
h,h '' 

Second, the value of the wave function at the origin is related to the leptonic decay width 
T{V — > e"'"e~) of the vector meson given by the following formula 



^ dz^l,v{z,r = 0) = where (0| Je^^iWI^^^) = ejv my e'' , (15) 



where fv is the coupling of the meson to the electromagnetic current and ey the isospin 
factor, which is the effective charge of the quarks in units of the elementary charge e: for 
the p meson it is the charge of the combination {uu — dd)/^/2, i.e. ey = l/\/2- Finally, one 
requires that the mean radius be consistent with the electromagnetic radius of the vector 
meson. 

In our calculation we consider the models proposed in Ref. [0] (which we refer to as 



DGKP) and |, 10| (which we denote NNZ(1994) and NNPZ(1997) respectively). 



Model DGKP The following wave function is adopted: 

4'\z, r) = 6,^^-, z{l - z) -^-^ f{z) exp(-u;VV2) • (16) 

V 47r V Nctv 



The parameter uj and the overall normalization are fixed by the condition ( [T^ ) and by the 
value of the leptonic width. (The exact values of these parameters as well as the form of the 
function f{z) are given in appendix A). Using the above form of the wave function ipy one 
can write p{r, Q) by replacing it in Eq. (^). In the case of model DGKP, it reads: 

p(r, Q) = 27rr2 J ^fy e, z(l - z)fiz)e-'^'^'/^ 2 z(l - z)QKo{er) . (17) 
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Models NNZ(1994) and NNPZ(1997) The wave function of the vector meson is given by 

^l;'^\z,r) = V-/^ V^ ^^^(i _ ^) K^(l - ^) - V.^ + ml]cP{r, z) . (18) 

The form of the radial wave function (j)(r, z) is given in appendix A. The expression for the 
overlap p{r, Q) then is 

p{r,Q) = 27Tr^^'^ J dz!^[m^yz{l - z) + ml] Ko{er) - eK,{er)-^^<P{r, z) . (19) 

The models NNZ(1994) and NNPZ(1997) differ only by the tuning of the parameters. 

Expressions ([T7| ) and (|l9| ) differ on several points. We study the influence of these differ- 
ences on the distribution p{r, Q) and on its integral N{Q) in the next sections. 



2.2 Comparison of the theoretical predictions with the HERA data at t = 

Before applying the models (|l6|,17) and (18,|l9) for the vector meson wave functions to our 
procedure of extracting the 5-matrix element, we additionally checked how they compared 
with the data on forward diffractive vector meson electroproduction at HERA. 

To this aim we compute da/dt\t=o using formulae (|l],§) and substitute the dipole-proton 
amplitude ^g^~^ by the one proposed in Golec-Biernat and Wiisthoff's saturation model^ 
{a{x,r) in the notations of Ref. |j|]). 

The theoretical curves corresponding to all three different models DGKP, NNZ(1994), 
NNPZ(1997) are plotted in Fig. ^ Let us stress that there is no additional tuning of any 
of the parameters of the wave functions. One of the models, NNZ(1994), is able to describe 
the data quite well, whereas the other ones underestimate the normalization by as much as 
a factor of 2. However, the energy dependence of the cross section is well described by all 
models. Qualitatively similar conclusions were reached in the study of Ref. fl^ . 

In order to understand this discrepancy between the different models, we plot in Fig. ^ the 
integrand p{r, Q)/r which appears in Eq. ( |lO|) together with the dipole cross section from the 
Golec-Biernat and Wiisthoff model. One can see that different models for the wave functions 
result in different shapes in r for p{r, Q)/r. It is clear from this plot that the exact shape of 
the function p(r, Q)/r for large values of r is very important in this case since it is weighted 
by the dipole cross section a{x, r) which is large in this regime. 



2.3 Overlap integral N{Q), mean radius rq and theoretical uncertainty on S{x,rQ,b) 

The results on the forward production cross section presented in Fig. ^ suggest that it is 
essential to study in detail the model dependence of the S'-matrix given by formulae (12) and 

dlD- 

Let us first study the properties of the quantity N{Q). In Fig. |^ we plot it for the three 
different models for the vector meson wave function. In the region of interest, 0.05 < < 
10.0 GeV^ they result in comparable N{Q) within 10 - 25%. This means that the average 

^ A similar study was done in detail recently in where different forms of the wave functions for po 
and J/'i/ were considered and used in the calculation of the forward longitudinal and transverse total cross 
sections. The meson wave functions were chosen according to general principles rather than just relying on 
existing models. 
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DGKP 


NNZ 


NNPZ 


[GeV^] 


rg [fm] 


A 


rg [fm] 


A 


rq [fm] 


A 


0.45 


0.35 


1.88 


0.49 


2.63 


0.52 


2.79 


3.5 


0.21 


2.24 


0.26 


2.77 


0.26 


2.77 


7 


0.16 


2.32 


0.20 


2.90 


0.19 


2.75 


27 


0.09 


2.49 


0.11 


3.04 


0.10 


2.76 



Table 1: The mean values rq of the dipole size for different obtained in models DGKP 
NNZ(1994) 1^ and NNPZ(1997) J7^. rq is evaluated according to formula (2C). Consistency 



with formula is checked by quoting the corresponding values of A, all quite close to 2 
and hardly dependent on Q^. 



^-matrix obtained using Eq.(12) is quite robust to a given choice of the model for the vec- 
tor meson wave function. In fact one can give some very general bounds on N{Q), nearly 
independent on the model for ipy From these upper bounds on N{Q) we could obtain an 
absolute upper bound for the S'-matrix which would be model-independent. We refer to ap- 
pendix B for technical details. We represent these bounds in Fig. ||. One can see that the 
models happen to be quite close to these limits (within say 25%). 

Second, we estimate the variation of rq with respect to the different models for the wave 
functions. We compute the mean value rq as: 

rq = e<'°s->- . (20) 

The results for the different models at different values of are given in Tab. 0. The values 
of rq are consistent within 25% for the various models. This confirms the validity of for- 



mula (13). 



The three models appear to be more or less equivalent for our purpose. For simplicity, 
we then only consider model DGKP [^ in the following. 

3 Impact parameter analysis of HERA data 

In this section we extract the 6-dependent 5-matrix from the HERA data. To this aim. 



we apply formulae (12,13) to the data, together with expressions (17,19) obtained from the 
discussion of the vector meson wave function in the previous section. We also try to deduce 
from these results the value of the saturation scale. 

3.1 Profile S{b) 

The experimental data for diffractive production of vector mesons are usually parametrized 
by the forward diffractive cross section da/dt\t=o and the logarithmic slope in momentum 
transfer B{t) as follows: 

dt dt \t^o 
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We take the available data |19, 2C] for the electroproduction of longitudinally polarized po 
vector mesons. These data are given for momentum transfer t below tmax = 0.6 GeV^, which 
enables us to determine reliably the S'-matrix only for impact parameter values larger than 



b = l/-v/tmax — 0.3 fm. In order to compute the Fourier transform appearing in Eq. (12), we 
have to assume an extrapolation of da/dt to larger values of t. The formula we use is 

r+oo \ 

+ / dtJo{hVi)^/£{f)\ . (22) 

Assuming that B{t) is non-increasing with t, as various sets of data seem to indicate, we 
choose different functional forms for the extrapolating functions S{t): 

• an exponential form exp( — d' t), with the constant d being of the order of B{txacix 

). This 

choice provides an upper bound on the 5-matrix for b close to 0. 

• a power law form t^". This choice is motivated by the fact that the data for photo- 
production at high t |21| indicate a i-dependence governed by such a power law, with 



a = 3. Such less steep t-dependence was also obtained in a theoretical calculation p4| . 
We consider here two values for a: a = 3 and a = 6. 

The parameters of these functions £{t) are fixed from the fit to the experimental points cor- 
responding to the highest values for i, in the region 0.4 GeV^ < t < 0.6 GeV^. These specific 
assumptions give an idea of the uncertainty on the determination of S{x,rQ,b) due to our 
lack of knowledge of the differential cross section for t > tmax- Another source of uncertainty 
is the experimental errors on the measurements of B(t) and da/dt\t=o- A complete error 
analysis would lie beyond the scope of this paper. Here we just estimate the influence of 
these uncertainties on S{x,rQ,b) by varying the measured quantities inside the 1-a error 
bars. By this method, we believe that we obtain a strict overestimate of the errors. The 
theoretical input N{Q) neeeded in formula ( p^ is computed within the model DGKP. The 
value rg is estimated using the procedure described in Sec. (|T|) and (|2|). The other models 
NNPZ(1997) and NNZ(1994) have also been tested and lead to very similiar results. 

We take three values of Q^: Q\ = 0.45 GeV^, Qj = 3.5 GeV^ and Ql = TGeV^. These 
values correspond to bins of the ZEUS analysis [|l9|, which we consider in the following. We 



note that the HI data |2C] lead to similar results. We recall that these values of correspond 
to the respective values of the dipole size rg (see Tab. ([l|)): rg^ = 0.35 fm, rg^^ = 0.21 fm, 
TQg = 0.16 fm. For these values of Q^, we take similar low values of x, namely xi = 4.7- 10~^, 
X2 = 4.3 • 10~^ and x^ = 5.8 • 10~^ respectively. The experimental slope B{t) is parametrized 
in the region t < 0.6 GeV^ by B{t) = Bo—ct, where the values of Bq and c are Bq = 9.5 GeV~^, 
c = 4.0GeV-^(forQ2 = 0.45 GeV^) and c = 6.1 GeV-^(for = 3.5 GeV^ and = TGeV^. 

The results for S{x, rg, b) are shown in Fig. (^, and are obtained using Eq. (|2^ ) applied to 
the ZEUS data. For each Q^, the curves corresponding to three choices for the extrapolation 
function £{t) are drawn. The uncertainty corresponding to the errors on the measurements 
themselves are computed for the intermediate value of Q^, and are represented in the lower 
part of the plot by a hashed band. 
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First, we observe that for b > 0.3 fm, the curves corresponding to different values of 
are clearly separated and ordered according to: 

S{xi,rQ^,b) < S{x2,rQ^,b) < S{x3,rQ^,b) . (23) 

As xi ~ X2 — 3:3, this means that the proton is less transparent to dipoles of larger size r. We 
also note that in this region of b > 0.3 fm, the results are not very dependent on the choice 
of extrapolation of the t-dependence of the data. However, for b < 0.3 fm (shaded region in 
Fig. (^)), the results become very sensitive to the asymptotic behaviour of the cross section 
da/dt, and data at larger t are needed to be able to make firm statements about the value 
of S{x, rq, b) near 5 = 0. We also note that the errors on the measurements for t < 0.6 GeV^ 
lead at most to an uncertainty of about 25 % for the S-matrix at impact parameter 6 = 0.3 fm. 



Finally, one can compute the total cross section for dipole-proton scattering using the 
following formula (which is nothing more than Eq. (P) averaged over r, for real S-matrix): 



(^tV(^> r))p = 2J d'h{l- {Six, r, b)),) . (24) 

Using formula ( [l2|) to replace {S{x,r,b))p and performing the integral over b, one sees 
that this cross section is determined by N{Q) and by the forward differential cross section 
da/dt\t=o- For photon virtualities Qf, Q2 and Q|, one obtains the values {cr1ot^)p ~ 14.9mb, 
10.6 mb and 7.5 mb respectively. These results are consistent with those of Ref. ||l], |l^, 11, 
12, 13, 14|, however, the relationship between and the radius rg quoted in these papers 
remains a large source of uncertainty in this comparison. 



3.2 Towards an estimate of the saturation scale 

One can try to estimate the quark saturation scale using the results presented in Fig. ^. 
From this analysis it is in principle possible to extract the dependence of this scale on the 
impact parameter b of the collision. We try the following phenomenological formula |23] for 
the S'-matrix: 

S{x,rQ,b) = exp{-rlQl{x,b)/A) . (25) 

This exponential form comes from a Glauber-like summation of multiple independent scat- 
terings of the quark-antiquark pair on the target nucleon |2^. It is also supported by 
the Golec-Biernat and Wiisthoff model |||] where the radius Roix) introduced there can be 
related to the saturation scale by Ro{x) = 1/Qs{x,b). However, in that work, it was assumed 
that the 5-profile of the S-matrix has the form of a sharp cutoff at a distance of order 1 fm. 



Using the approximation rq ~ + rriy one obtains ~ 1 — 1.5 GeV^ and Qg ~ 

0.2 GeV^ for b = 0.3 fm and b = 1.0 fm respectively. These values are consistent in order 



of magnitude when computed at different Q^, which confirms the relevance of formula (25). 
These results would suggest that at small impact parameter values the saturation scale is in 
the semi-hard regime and support the onset of the shadowing effects at HERA. One should 
however stress that these estimates are rather rough and should be taken with care since they 
strongly depend on the precise value of rq. Indeed, rq comes squared in the formula for Q^. 
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4 Summary 



In this paper we have shown that using the HERA data on diffractive electroproduction 
of vector mesons it is possible to extract the S'-matrix element S{x,r,h) for dipole-proton 
scattering in impact parameter space. By considering the full t dependence of this process 
we have shown how to obtain the ^-matrix element averaged over dipole sizes r. Due to the 
particular properties of the overlap function between photon and meson wave functions this 



average can be replaced by the S-matrix element evaluated at r = tq 2/ y Q'^+rriy. We 
have then shown that our results on the S-matrix are only weakly dependent on the choice for 
the model of the meson wave function. Since the data are available only for low momentum 
transfer, t < 0.6 GeV^, we have shown that our results are reliable for impact parameter 
values larger than b ~ 0.3 fm. It would be very interesting to extend the measurement of 
diffractive electroproduction of vector mesons to larger values of t. This would enable us to 
explore the very interesting regime of central impact parameter collisions at HERA. We have 
also estimated the value of the dipole cross section integrated over the impact parameter and 
have found it to be consistent with other analyses. Finally we have discussed how to use this 
result to estimate the saturation scale at HERA collider. We have suggested that the onset 
of the shadowing effects can be dependent on the impact parameter of the collision and that 
for central collisions it could occur in the semihard regime. However further theoretical and 
experimental studies are necessary. 
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A Models for the meson wave function 



In this appendix, we detail the phenomenological forms we adopt for the meson wave function. 

A.l Model by Dosch, Gousset, Kulzinger, Pirner (DGKP, Ref. [§]) 

The form of the vector meson wave function is given by: 



i^'v\z,r) = 5,^_-,^Au:{z{l-z)f'^ exp {-\'^{z-l/2f\ exp {-\u:\^ ] . (26) 



The parameters chosen are the following: 

Lo = 0.330 GeV and N = 4.48 . (27) 

An additional feature of this model is that the mass niq of the quarks that appears in the 
photon wave function is running with This gives important effects only for very low Q 
and for photoproduction. In this regime, it is argued that the quarks should have constituent 
mass. The formula for the running mass reads: 

r0.220GeV-(l-gVQo) for < q2 ^ ^ 05 GeV^ 
y for > Qo 

This model is referred to as DGKP. 

A.2 Models by Nemchik, Nikolaev, Predazzi, Zakharov (NNZ(1994), Ref. ^ and NNPZ 
(1997), Ref. fligi) 



In this model, the radial wave function (t>{r,z) appearing in eq.(|l9|) satisfies the following 
normalization condition: 



1 = § ^2^^ J d'rl^rnyiz, r) + {z' + (1 - zf)idMz, r))^} . 



(29) 



This equation is nothing else but the normalization condition (U) , applied to the transversely 
polarized meson wave function. The obtained normalization factor ^'o(15') is assumed to be 
the same for a longitudinally polarized meson. The function (j){r, z) is defined by: 

^{z,r) = ^o{lS) 



} 



A{z,r)B{z,rf 

The functions A and B are given by 



+ ^^^y ^. TKi{r.A{z,r)/B{z,T)) \ . (30) 



A\z, r) = 1 + / \ - AC^a\r)mi 

z{l — z) ^ 



d2, ^ C^a^r) 



(31) 



13 



and a(r) = 3/(8mgQ<j(r)). The following prescription is taken for the running coupling as(r): 

47r 

as{r) = ao for r > and as{r) = — — „ „ for r < , (32) 

/3olog(l/AV^) 

where the parameter values are 

rs = 0.42 fm, oq = 0.8, A = 200 MeV . (33) 

The masses of the quarks are taken to be 0.15 GeV. 

The other parameters ^'o(15'), C and are chosen taking into account several con- 
straints: the normalization condition (|29| ) for the wave function has to be satisfied, and the 
value of the leptonic decay width must agree with the experimental measurement. Addition- 
ally, the mean radius of the meson has to be of the order of a hadronic scale. Two different 
sets of parameters are considered: 



C = 0.25, = 0.76 fm^ (reference @) 
and C = 0.36, = 1.37 fm^ (reference |l 



(34) 

These two choices are refered to as NNZ(1994) and NNPZ(1997)respectively. 
B Model-independent bounds on the overlap integral N{Q) 

In this appendix we explore more the properties of N(Q) and show that upper bounds can 
be given on this quantity, regardless the model adopted for ipy- These provide upper bounds 
on S{x,rQ,b). 

N{Q) can be seen as a scalar product of the two wave functions, N{Q) = (ip^lipv), see 
Eqs. (^) and (10). The Cauchy- Schwartz inequality then applies. It leads to an upper 



estimate for the integral: the product of the two wave functions is smaller than the square 
root of the product of the integrals of the squared wave functions. Using the normalization 
condition ( [l^ ) for the meson wave function^ and computing the full integral of the squared 
photon wave function, one eventually obtains: 

^,(g)<,,,^[l_«^,,49|-^.a„,r.^| . (35) 



The r.h.s. of eq.(|35|) grows with and thus this inequality is only interesting in the small-Q^ 
region. Although independent of the precise form of the vector meson wave function, this 
bound nevertheless depends on the masses of the quarks. 

Second, with the additional assumption that ipviz, r) is maximum for r = 0, we can write: 
N{Q)<2ev j d^rdzV7*(-z = l/2,r;Q)^y(z,r = 0) . (36) 



* Strictly speaking, this bound is only valid for models for which the condition is enforced for the 
longitudinally polarized vector meson. 
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The integrals over z and r in the r.h.s. are now factorized. The integration over r is performed 
analytically, while the one over involving ^/^y , can be expressed as a function of the coupling 
fy of the meson to the electromagnetic current, using the relation (p^). This finally leads to: 

iV(Q) < 2V^/y ^2 + 4^2 • (37) 

The r.h.s. vanishes like 1/Q at large Q, which makes this inequality useful for large Q. 

The region in the {Q'^,N{Q)) plane which is forbidden by these bounds is depicted in 
Fig. (|5|). One sees that all the models are very close to the upper bounds. 
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Figure 1: Diagrammatic representation of diffractive production of vector mesons. 7*(^) is 
the virtual photon of momentum q, V{q + A) is the produced vector meson with momentum 
q + A. The target proton is scattered elastically. 
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r(fm) 



Figure 2: Logarithmic distribution of dipole sizes at the photon vertex p{r,Q)/N(Q), for 
different Q^. 
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Figure 3: The forward longitudinal cross section for the production of po vector meson. Data 
are from |jl^ and the three curves correspond to different assumptions for the wave functions 
of the vector meson. Solid curve: model DGKP Q, dashed curve: model NNPZ(1997) ||l^ 
and dotted curve model NNZ(1994) ||9|. The dipole cross section is taken from Golec-Biernat 
and Wiisthoff model 



18 




Figure 4: The overlap of the wave functions of the photon and vector meson p(r, Q) plotted 
as a function of r. Solid curve: DGKP Q, dashed curve NNPZ(1997) ||l^, dotted curve 
NNZ(1994) g. The dashed-dotted curve corresponds to the dipole cross section by Golec- 
Biernat and Wiisthoff [|^ calculated for x = 10~^. The relative normalization of the dipole 
cross section and the wave functions is not conserved. 
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Figure 5: Integrated dipole distribution N{Q) as a function of Q^. The curves correspond- 
ing to the different models for the meson wave function are shown. The shaded area is 
theoretically forbidden. 
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Figure 6: S'-matrix for dipole-proton scattering as a function of impact parameter b. Three 
different are considered, corresponding to three typical values for the size of the interacting 
dipole (which are estimated according to rg = e*-'°s'^>p). For each value of Q^, the curves 
corresponding to three extrapolations of the data for t > 0.6 GeV^ are shown. The shaded 
band indicates the region of impact parameter b where the choice of this extrapolation is 
crucial, and thus where our extraction is not reliable. The hashed band on the bottom is an 
estimate of the errors due to the experimental uncertainties on d(j/dt\t=Q and on B{t) (for 
t < 0.6 GeV^). 
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